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Abstract 

A master equation describing the evolution of averaged molecular state occupan- 
cies in molecular systems where alternation of molecular energy levels is caused 



by discrete dichotomous and trichotomous stochastic fields, is derived. This study 
is focused on the kinetics of quasi-isoenergetic transition processes in the pres- 
ence of moderately high frequency stochastic field. A novel physical mechanism 
for temperature-independent transitions in fiexible molecular systems is proposed. 
This mechanism becomes effective when the conformation transitions between quasi- 
isoenergetic molecular states take place. At room temperatures, stochastic broad- 
ening of molecular energy levels predominates the energy of low frequency vibra- 
tions accompanying the transition. This leads to a cancellation of the temperature 
dependence in the stochastically averaged rate constants. As examples, physical 
interpretations of the temperature-independent onset of P2X3 receptor desensitiza- 
tion in neuronal membranes, as well as degradation of PER2 protein in embrionic 
fibroblasts, are provided. 
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1 Introduction 



Transient processes in physical, chemical, and biological systems reflect the 
evolution of state occupancies associated with different types of electronic, 
nuclear, and spin degrees of freedom. From the physics point of view, the de- 
velopment of the microscopic models that are able to shed light on possible 
mechanisms of transient processes, is an important task. Biological systems 
are characterized by complex molecular structure with fast, intermediate and 
slow motions of separate molecular groups. Most transient processes in biosys- 
tems can be studied in the framework of phenomenological approach that is 
based on the equations with certain set of kinetic parameters (rate constants). 
In these cases, however, any experimentally observable transition process oc- 
curs simultaneously with the more fast transitions. For instance, the charac- 
teristic times of conformation transitions in biopolymers (nucleic acids and 
protein chains) are several orders greater than the characteristic times of the 
vibration relaxation. Therefore, the phenomenological parameters including 
the occupancies of conformation states and the rates of transitions between 
them, appear averaged over the much faster motions. Numerous experimental 
studies show that in biosystems, the temperature dependence of the reaction 
rates exhibits the Arrhenius like behavior, i.e. an exponential increase of the 
rate with temperature. However, there are cases when a temperature depen- 
dence smoothes away or becomes anomalous (opposite to "normal"). Several 
attempts were made to explain these effects in the framework of phenomeno- 
logical approach, without exploring the microscopic properties of molecular 
structures [TpfSPfS] . The goal of this study is to provide a quantitative mi- 
croscopic description of both the smoothed away (temperature-independent) 
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and normal (temperature-dependent) cases for the multi-type reaction path- 
ways arising within the complexity of molecular structure. 

In biology, there are three known types of reactions where microscopic effects 
could be essential. The first type is associated with the reduction-oxidation 
processes caused by the photoinduced electron transfer in the photosynthetic 
units of plants and bacteria. The basic mechanism of these reactions involves 
electron tunneling between rigid redox centers like metal-containing porphyrin 
rings and quinones. The tunneling mechanism is supported by the absence of 
temperature dependence for the transfer rates in wide temperature regions 
including the room temperatures |6f7f8] . The physical model describing an 
electron transfer between redox centers is well known, and it is associated with 
the existence of strong coupling with local vibrational modes. The respective 
analytic expressions have been derived by Jortner [QfTO] who has general- 
ized the classic Marcus's theory [TT] of transfer processes for a quantum case. 
When contact conformations between the centers can be stabilized, an elec- 
tron recombination rate averaged over the conformation states, has anomalous 
temperature dependence demonstrating a certain rise with the temperature 
decrease [T2] . 

The second type of reactions is associated with the biochemical catalysis in 
the specific protein or metal-protein centers. As a rule, such catalysis leads 
to a notable structure reorganization in the active centers as well as breaking 
of the old and creation of the new chemical or hydrogen bonds between the 
respective structure units. Catalytic process consists of various chemical steps 
that include not only an electron and a proton transfer, but a two-electron, 
an electron-proton, a hydride or a hydrogen transfer as well [T3|Ti] . Conven- 
tionally, all such reactions are described in the framework of enzyme kinetics 
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|15] where the rate constants are calculated using the concept of an activated 
complex strongly connected with the corresponding reaction coordinate. The 
motion of the reactants along reaction coordinates is accompanied by over- 
coming the activation barriers, so that the rate constants experience an acti- 
vation temperature dependence similar to the Arrhenius's one. However, in the 
framework of Marcus's theory, a barrierless behavior becomes also possible, if 
only the value of the reorganization energy would coincide with the driving 
force of the reaction (due, for instance, to a pronounced nuclear displacement 
in rigid complexes), thus cancelling the exponential temperature dependence. 
Such situation is probably realized during the binding of the NO molecule to 
the protoheme. The binding appears as a kinetically mono-exponential pro- 
cess with approximately temperature-independent rate constant in the region 
of 200-290 K [16]. 

The third type of reactions is associated with the conformation transforma- 
tions. Some of these reactions do not exhibit the Arrhenius's temperature 
dependence even in the range of room temperatures. For instance, it has 
been demonstrated that the closing rate for both the single-stranded DNA 
and the RNA hairpin-loop fluctuations, as well as the rate for the cyclic /3 
-hairpin peptide folding, are weakly dependent on temperature in the inter- 
val from 10 to 60°C showing zero or even slightly negative enthalpies pT|ITS] . 
Another example of the temperature-independence in the gating kinetics of 
membrane proteins has been recently provided for the desensitization onset 
of P2X3 purinoreceptors [I9]. The kinetics of this process is characterized by 
the two rate constants, each being temperature-independent in the range of 
25-40°C. A very recent example is the observation of a single-exponential ki- 
netics of temperature-insensitive (in the range of 27-37°C), period-determining 
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process in the mammalian circadian clock [20]. Note that, in these cases, the 
experiments cover a physiologically important temperature range. 

While physical mechanisms explaining both the temperature-independent elec- 
tron tunneling and the barrierless ligand binding are more or less clear for 
two noted types of reactions, and can be well recognized in the framework 
of Marcus- Jortner approach, it is not the case for the third one. Indeed, the 
conformational transformations in flexible molecular structures should involve 
activation of low-frequency vibrations accompanying the quasi-isoenergetic 
configuration transitions. This would make the temperature dependence of 
the transition rates exponential. However, several experimental results have 
demonstrated a presence of the temperature-independent reaction pathways. 
Therefore, to describe the quasi-isoenergetic transient processes, the theoret- 
ical approach has to be modified. 

In general, the current theoretical description of transient processes in flexible 
molecules as well as in molecular chain structures employs two approaches. 
One is based on quantum mechanics and uses the recent ab initio methods 
(mainly, a density functional theory [H]) along with the ab initio molecu- 
lar dynamics simulations [22]. Another approach involves the development of 
various physical models that allow one to describe the amplitude and the ki- 
netic characteristics of the transient processes in terms of definite physical 
parameters like vibration coupling, energy bias, fluctuation frequency, level 
broadening etc. [23j. Most of these parameters already appear in molecular 
Hamiltonians, others are described with the rate constants of the molecular 
transitions (see e.g. [2^|25] and refs. there in). The problem is, however, that 
dealing with large molecular systems like proteins or nucleic acids in aqueous or 
specific biological environment using ab initio methods is very expensive com- 
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putationally. For instance, it takes 4 days using a cluster of four Linux-based 
computers to calculate a nanosecond molecular dynamic trajectory for a small 
protein containing ~5000 atoms and surrounded by ~15000 water molecules 
[26]. An overwhelming majority of biomolecular transient processes requires 
times ranging from milliseconds to seconds. So, correct computer simulation 
of these processes is currently inaccessible, even by modern supercomputers, 
even in times from thousands to millions of years. 

Meanwhile, the rigorous approach of statistical mechanics allows to describe 
the evolution of an open (interacting with the environment) molecular system 
in the formalism of nonequilibrium density matrix. This formalism has been 
successfully applied to describe exciton and charge transfer in different types 
of molecular systems P7|28f29] . The fundamental advantage of the density 
matrix formalism lies in the unified description of transfer processes by taking 
into account both the heat bath vibrations and the stochastic energy fluctua- 
tions [5U|31|32] . Furthermore, contribution from the stochastic effects can be 
analyzed exactly in a non-perturbation manner [30|32f33f34f35j . Another ad- 
vantage of the noted formalism is that it allows one to reduce the description 
of evolution behavior to the finding of the solution of a set of kinetic equations 
for state occupancies with precise rate constants. Therefore, the interpretation 
of the results in terms of the rate constants becomes clear and useful for both 
physicists and biologists. 

In this paper, we develop an approach that is specially adapted for descrip- 
tion of transitions between quasi-isoenergetic molecular states, with use of 
nonequilibrium density matrix method. We argue that, for quasi-isoenergetic 
transitions, both low-frequency molecular vibrations and stochastic alterna- 
tion of molecular energies by thermodynamic fluctuations, should be simulta- 
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neously accounted for. We obtain rather uncomplicated expressions for tran- 
sition rate constants in various limiting cases. We demonstrate that in these 
cases, different types of temperature behavior - from exponential tempera- 
ture dependence of Arrhenius like type to the temperature independence, are 
adequately reproduced. The examples include a physical explanation for the 
temperature-independent onset of P2X3 receptor desensitization in neuronal 
membranes [19] and temperature-insensitive, period-determining process in 
mammalian circadian clock [20] . 

2 Model and basic stochastic equation 

In molecular systems, any transition process contains contributions from differ- 
ent types of vibrations. Some of them define the common electron- vibrational 
molecular states involved in the transitions while the others can be referred 
as a heat bath. Just owing to coupling with the heat bath, the transitions be- 
tween the molecular states appear in the form of kinetic process. This process 
reflects an evolution of the molecular system towards the equilibrium, where 
state occupancies satisfy the Boltzmann's relations. As it is well known from 
the physics of transfer processes in condensed matter, the stochastic fields 
created by interior motions in a dynamic system, can play an important role 
in evolution of state occupancies. Direct introduction of these fields into the 
Hamiltonian of dynamic system significantly simplifies the consideration of the 
complex motions in the system. This allows one to use only a limited number 
of general characteristics like amplitudes and frequencies of the stochastic field 
(see, for instance, refs. [32|33f34f35f36] ). The stochastic approach clarifies the 
physics of transition processes on the semi-phenomenological level. It becomes 
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especially important in the case of biological macromolecules where a number 
of different types of nuclear motions, including the stochastic motions, are 
observed. It has been experimentally shown [37J that structural movements of 
separate polar residues in proteins create high frequency vibrations, so that 
specific electric field fluctuations are present in a large number of biologi- 
cal macromolecules. Direct observation of fluctuation movements in biologi- 
cal molecules has been done with use single-molecule spectroscopy. Namely, 
the observation that the fluorescence of a single flavin adenine dinucleotide 
cofactor at room temperature turns on and off indicating the stochastic en- 
zymatic turnovers of a flavin reductase [38|l39] . The energy fluctuations at 
ambient temperatures have been also observed in a photosynthetic peripheral 
light-harvesting molecular complex of type 2 00111] ■ Thus, there exists an ex- 
perimental evidence supporting a stochastical approach to describe transition 
processes in biological macromolecules. The simplest way for such descrip- 
tion is to incorporate the random processes in a stochastic Hamiltonian and, 
therefore, in the stochastic Liouville equation. 

2. 1 Stochastic molecular Hamiltonian 

When the exact spatial molecular structure is unknown or uncertain (often in 
the case of biological macromolecules), the modeling of the transition process 
has to be based on general physical principles. While developing such model, 
two different types of motions in a macromolecular system have to be consid- 
ered. First, vibrations in a heat bath. Due to phonon exchange with a heat 
bath, the transitions between the molecular states appear as the relaxation 
process where the phonon energy covers the difference between the energies of 
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these states. Second, high frequency vibrations. The energy of the correspond- 
ing phonons exceeds the noted energy difference signigficantly. This makes a 
phonon exchange ineffective. However, high frequency motions create inte- 
rior random fields that alternate the molecular energies. The above examples 
show that at such alternation, the molecular energies can be stochastically 
time-dependent and, thus, the transitions between the molecular states (ac- 
companied by creation and annihilation of low frequency phonons) occur on 
the background of randomly alternated energy levels. Taking this notion into 
account, we consider a model where molecular energy levels exhibit stochas- 
tic alternations caused by interaction with the surrounding structural groups. 
Rather, the transitions between the energy levels are associated with the har- 
monic vibrations (phonons) of a heat bath. Respective Hamiltonian of the 
entire system (molecule -|- bath) can be represented in the form 

HenUreit) = E {^"(t) + E Mf^X + /^a) + hw^i^t + 1/2)]} \a) {a\ 
a X 

+ J2yab{l-Sa,b)\a){b\ (1) 
a,b 

where Ea{t) — E^^^+AEa{t) is the energy of the ath molecular state {E^^^ is the 
energy in the absence of an external time-dependent field, AEa{t) is the energy 
variation caused by the stochastic fields), Kax is the parameter that charac- 
terizes the deviation of nuclei along Ath normal coordinate from the nuclear 
equilibrium position, and (jJx is the Ath vibration mode (phonon) belonging to 
the bath Hamiltonian. Quantities ^^I^ and P\ refer to the operators of creation 
and annihilation of a phonon in a heat bath. Thus, the entire system under 
consideration consists of two main parts. The first one is the dynamic quan- 
tum system (S) with the set of molecular states {\a)} and respective energies 
{Ea{t)}. The second part represents a heat bath (B) with set of vibrational 
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modes {u)x}. The transitions between the molecular states are determined by 
the term containing the matrix elements Vba- In what follows, we represent the 
main parts of the Hamiltonian ([T]) in completely diagonal form. To this end, 
it is convenient to use the Holstein's polaron-like transformation |12|^ while 
defining the respective unitary matrix in form U = exp (^J2a chalet) with 
f^a = Y.X di^^ {(^x — l^x) being the operator of the nuclear displacements in the 
ath molecular state {g^''^ = Kax/hjJx is the dimensionless coupling). Multiply- 
ing the Eq. ([1]) from the left by U and from the right by we transform the 
starting Hamiltonian ([T]) to 

H{t) = Hs{t) + HB + V . (2) 

Here, 

Hs{t)=Y.{E^{t) + ^E'f''^^)\a){a\ (3) 

a 

is the molecular Hamiltonian where a proper energy of the ath state in- 
cludes polaron shift AE^^"'-^ = —Y.x '^oa/^a- Bath Hamiltonian coincides 
with Hamiltonian of harmonic vibrations (phonons) , 

Hb = Y: huJx{PtPx + 1/2) • (4) 

A 

Off-diagonal part of entire Hamiltonian appears now in the form 

\/ = ^(l-MKb|a)(&| (5) 

a,b 

where Kb = Kfec'^'"'' and aab = c^a - (^b = Y.x g^ab {Px - /^a)- Operators Vab 
characterize the transitions between different molecular states when the tran- 
sitions are accompanied by the nuclear displacements (respective couplings 
are g^^^^ = g[^^ — gh^^)- Thus, the transitions between the molecular states are 
performed with the " phonon-dressed" off-diagonal interactions. 
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2.2 Stochastic equation for a diagonal part of density matrix 

To perform the analysis of the experimental data, a master equation for de- 
scription of evolution of the observable state occupancies Pa{t), is reuired. 
Since the molecular energy levels Ea{t) are the stochastic quantities, this mas- 
ter equation appears as the equation averaged over random realizations of 
the energy shifts AEa{t). We denote the averaging via the symbol ((...)) so 
that Pa(t) = {{Vait))). Here, Va(t) is the nonaveraged state occupancy which 
is defined as Va(t) = {a\p(t)\a) with p(t) = trs Ps+sif) being the molecular 
density matrix. This matrix is determined as a trace over the bath states on 
density matrix ps+B{t) of the entire quantum system including the heat bath. 
The density matrix method assumes that irreversible (relaxation) transitions 
between the quantum states of an open dynamic quantum system S 
sociated with an energy exchange between the system and the heat bath (via 
creation and annihilation of phonons). Since the coupling between the molec- 
ular states and the heat bath is concentrated in off-diagonal interaction ([5]), it 
becomes important to derive a master equation for molecular occupancies in 
the form that allows one to control an expansion over the noted off-diagonal 
interaction. Due to the exact relation {a\p{t)\a) = {a\pd{t)\a) , an equation for 
diagonal part of molecular density matrix, pd(t), is quite sufficient. The de- 
tails of the derivation procedure can be found, for instance, in refs. pHl25P^ . 
Following this procedure while taking into account that, in the considered 
case, the molecular energies are stochastically time- dependent quantities, the 
master equation for the Pd{t) should appear in the stochastic form. 

The starting point is the stochastic Liouville equation 

PsMt) = -imPsMt) (6) 
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where £(t) = (l/h) [H(t),...] is the Liouville superoperator associated with 
the Hamiltonian ([2]). In hne with the structure of this Hamiltonian, the C{t) 
consists of two Liouville superoperators, Co(t) = (l/h) [HQ(t), ...] with HQ(t) = 
Hs{t)+HB and Cy = (1/^) [V, ...] so that = CQ(t) + Cv Let introduce the 
projection superoperators and Tnd = I—Td- These will expand any operator 
into its diagonal component (the matrix elements are the occupancies) and its 
off-diagonal component (the matrix elements are the coherences). The action 
of the operators and Tnd on the Liouville equation ([6]) generates a coupled 
set of differential equations for diagonal and off-diagonal parts of the density 
matrix, p^s+si^) = TdPs+B{t) and pf+sit) = fndPs+B{t), respectively. These 
equations read 

pf%{t) = -ifnd mp^sikt) - ^cyp^sUt) . (7) 

After substitution of the second equation into the first one, the following 
integro-differential equation for the diagonal part of the density matrix is 
achieved 

P^sUt) = -f dt' [fd CyS{t, t')Cyp's'Ut')] (8) 


where 

S(t,t') = e~'-^t''^^^"^^(-^ (9) 

is the stochastic evolution superoperator. Bearing in mind the definition of 
the molecular density matrix, one arrives to the equation 

Pd{t) = - dt' trB[fdCyS{t,t')Cypf_,B{t')] ■ (10) 

In many practical cases, it is common that a characteristic time of transition 
processes, Ttr, significantly exceeds the characteristic time of the vibration re- 
laxation, Tmb- This makes it possible to refer the vibrations to the heat bath. 
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Moreover, the validity of the inequahty Tyn, <^ Ttr allows one to use a factoriza- 
tion p^'slsit) = Pd{t)pB where ps = exp{-HB/kBT)/trB exp{-HB/kBT) is 
the bath equilibrium density matrix {ks and T are the Boltzmann's constant 
and the absolute temperature, respectively). Thus, if the transition processes 
occur on the time scale of At ^ r^b, one can describe these processes with 
the coarse-grained master equation 

Pd{t) = - fdt'M{t,t')p,{t') (11) 
Jo 

where 

M{t, t') = tTB [fd CySit, t')CvPB] (12) 

is the stochastic transition superoperator. 



3 Master equation for the averaged molecular occupancies 

Eq. ffTTl) can be treated as the basic coarse-grained stochastic equation for 
the rigorous description of the transitions between the molecular states on the 
time scale At ^ T^ib. Its form is quite suitable for an expansion procedure. 
Thus, if a description is restricted by Born approximation over the V one has 
to set V = in the evolution superoperator Qj. It yields 

Pd{t) = -{l/h^) J^t' trB{f4V,U{t,t')[V, pd{t')pB]U+{t,t')]} (13) 

where we have introduced the zeroth order evolution operator U{t, t') = 
Texp [—{i/h) jl, dT{Hs{T) + Hb)] {T is the Dayson's chronological operator). 
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3. 1 Stochastic kinetic equation for the state occupancies 



In what follows, we concentrate on the transition processes caused by a weak 
coupling of the molecular states to the bath vibrations (phonons) so that the 
" dressed" transition operators Vab in off-diagonal interaction (jS]) are considered 
as perturbation (see also refs. [21|32f45|H5] ). Using the definition V{a;t) = 
{a\pd{t)\a), one derives the following stochastic equation for the occupancies 



V{a;t) = -Y. / dt'\gabit,t')Via;t')-gUt,t')Vib;t'] 



(14) 



where kernel 



Gab{t, t') 



VbnrRe 



e*''-(*~*')7^„,(t-t')-^afe(t,0 



(15) 



exhibits a stochastic behavior through the stochastic functional J-'ab(t,t') = 
exp {(t/h) ll, dT[^Ea{T) - ^Eb{T)]}. In Eq. ^ab = (l/7i)[(i?f + i?^'))- 
{E^'^ +E'f^'^y\ is the transition frequency. Coupling to the heat bath is reflected 
in the correlation function IZabiT) = trs (^Pb exp [crafe(0)]exp [crab(r)]^. If the 
bath is associated with the harmonic oscillations, then cxabir) = J2x 9ab [bt^^^'^^'^- 
bxe''^^''] and, thus 

n,b{r) = e-^'^'W (16) 



where quantity 



coth 



'1 — coslot) + i sincjr 



(17) 



is of basic importance for any calculations done on the spin-boson model 
[29]l3MITrE8B9]l5^ In Eq. ([E 
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is the spectral function that includes an information on both the vibrational 
structure of macromolecule and the characteristics of coupling to the bath 
vibrational modes. 

Form (fT7|) is quite suitable for estimation of rate constants in the case of a 
strong coupling to the bath. This work, however, is focused on the transi- 
tion processes where a weak coupling is realized between the molecular states 
and the phonons. Therefore, it is more convenient to submit the correlation 
function in the form 

nabir) = e-'''^'' Qab{r) (19) 

where quantity Dab = (S'ifc'')^ {'^n{uj\) + 1) specifies the Debye- Waller factor 
while function 

Qabir)= rfcue-n E ^Ip.i(^a) VtXt) '^[^-Epac^a] (20) 

fixes the time dependence of the correlation function. In (l20l) . Ip{z) is the 
modified Bessel function, and n{uj) = [exp (hu/kBT) — 1]^^ is the Bose dis- 
tribution function. Coupling to the phonons is concentrated in the param- 
eter zx = 2{gl^^)'^ ^n{ijjx)[n{ujx) + 1], while index px indicates the number of 
phonons of the Ath mode that accompany the transition. The minimal number 
of phonons is equal to 1, so that J2\ P\> 1- 

Supposing the small nuclear displacement along the Ath normal coordinate, 
we consider a behavior of the correlation function at ^ 1- This allows one 
to employ the asymptotic form I^p^^(zx) ~ (2:^/2)1^^1 /IpaI') which indicates that 
the main contribution to transitions comes from the single-phonon processes. 
Therefore, by setting /o(z ^ 1) ~ 1 and Ii{z 1) ^ {z/2), one can see that 
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the form (120!) reduces to 

Qa.(r) = E(^?i?)'^A(r), (21) 

A 

where Rx{t) = n{oJx) exp {iuj\r) + [n{uj\) + 1] exp {—iuj\T) is the single-phonon 
factor. Bearing in mind that, at weak nuclear displacements, the role of the 
Debye- Waller factor is not important (exp {—Dab) ~ 1)? one can further specify 
the stochastic equation (fT4|) setting IZabit — t') ~ Qab{t — t'). 

3.2 Kinetic equation for the averaged occupancies 

To derive the equations for the observable occupancies Pa{t), one has to av- 
erage the stochastic equation fll4p . Here, we consider only the case of fast 
variations for each random energy difference AEair) — AEbir) , when the char- 
acteristic times of the stochastic and the transition processes {rstoch and Ttr, 
respectively) obey the inequality Tgtoch Ttr- This inequality allows us to use 
a decoupling procedure {{J'ab(t,t')Va(t'))) = {{J^ab{t,t')))Pa(t'). Moreover, a 
non-Markovian character of the transition process becomes not important. Ac- 
tually, let us set t = t' -|- r and rewrite Pa(t') as Pa(t — r) = exp {—T-^)Pa(t). 
On the time scale At ~ Tstoch the factor exp (— t^) can be estimated as 
exp {—Tstoch/Ttr) ~ 1 and thus Pa{t') ~ Paif)- This is also supported by the 
conclusion that in the Born approximation over the V [51] , the delay processes 
do not introduce any noticeable contribution to the transition kinetics [52] . 

The estimation of the averages like {{J^ab{t,t — r))) = Fabir) is given in nu- 
merous papers (see, e.g. refs. [3T|32ll33f36f53j ). The results indicate that inde- 
pendently of a specific form of quantities Fab{T), all of them decay with their 
own characteristic times, Tstoch- Taking this fact into account and using the 
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condition Tstoch ^ Ttr, one can shift the upper hmit in the integral of non- 
Marcovian stochastic equation (!T4l) from t up to infinity. Thus, providing the 
averaging over the fast random reahzations, we transform this equation into 
a balance like kinetic equation, 



To derive an analytic form for the rate constant, one has to calculate function 
Fab{T). To this end, let us consider the Kubo- Anderson random process [5l|55] . 
During such process, the jumps from one value to other happen at random 
time and independently from one another with the average frequency u. In our 
case, a stochastic alternation is associated with energy differences Ai?a(r) — 
AEb^r) = haahij). We limit ourselves by the case when fluctuations of the 
aab{T) are independent of the precise molecular states a and h. This situation 
is adequate for transitions in two-level systems and in flexible molecular chains 
containing identical structure groups. If aabir) = a(r), then Fabir) = F{t) = 
((X(r))) where X{t) = exp [i Jq dr'a^r')] is the stochastic functional. Below, 
two examples of Kubo- Anderson process are given. 

Dichotomous process (DP). During the DP a random quantity fluctuates be- 
tween two equiprobable values +a£, and — o"d with a mean frequency V£, so that 
oiab{t) = otDp{t) = ±o"_D- Bearing in mind the exact properties of the DP, so as 
«z)pW = o^D) ((«Dp(^))) = 0) and {{aDp{t)aDp{0))) = exp (-z/z)|t|), and 
using the Shapiro-Loginov theorem [36] in the form of relation d{{aDp{T)X (t))) / dr = 
ih'Dd{{X (r))) / dr + ia'j){{X (t))) , one obtains the following exaci equation (see 



(22) 



b 



with the averaged transition rate constants 




(23) 
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F(r) + udF{t) + alF{T) = . 



(24) 



This equation is similar to the one for a harmonic oscillator with frequency 
an and damping parameter Pr). Respective solution reads 

F(r) = (A:ie-^^" - k2e-'''^) / {h ~ . (25) 

The quantities 

= 1^0/2 ± v^(z/d/2)2 - al (26) 

determine the above noted characteristic stochastic times t^I^^\ = (Refci,2)~^- 
At high fluctuation frequency, owing to condition z/|) ^ Aajj the rate ki 
strongly exceeds the rate ^2- Therefore, at the time scale At ^ TstLh^ 
characteristic function appears as F{t) ^ exp (— 7£)r), where 

7d = (Tl/un (27) 

exhibits itself as the only damping parameter specifying the characteristic 
stochastic time r^toch = 1d^- An approximation of random energy variations 
by the DP is well justified for a two-level system. 

Trichotomous process (TP). The random energy difference AEa{t) — AEh{t) = 
haab{t) exhibits itself in the peculiar form of Kangaroo process [361153] where 
the aabit) fluctuates between the three equiprobable values, aab(t) = aTpif) = 
0, ±(Jt and {{aTp{t))) = 0. At the TP, the fluctuation frequency ut is assumed 
to be independent of the values of the amplitudes and, thus, the TP represents 
a member of the Kubo- Anderson processes. In this case, the Shapiro-Loginov 
theorem for the rules of differentiation of the kth order correlations reads 
rf((4p(t)X(t)))M = -z.T((4p(t)X(t)))+z/T((4p(t)))((X(t))) + ((4p(t)X(t))). 
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Using this theorem and bearing in mind the exact properties of the TP so 
as {{aTp(t)aTp{0))) = (2(7^/3) exp(—z/T|i|) and «|^p(t) = cr^aTpit), for the 
characteristic function F{t) one obtains the following exact equation 

F{t) + 2utF{t) + (4 + al)F{T) + (2z/T4/3)F(r) = . (28) 

This equation indicates the presence of three characteristic stochastic times. 
But, if the condition (z/|. + cr^)'^/^ ^ z/Tcrf. is valid, then on the time scale 
At ^ (z/|, + (t|.)^^/^ the characteristic function manifests itself as a single 
exponential drop so that F{t) ^ exp (— 7T'r). The decay parameter 

It = ^l/vT (29) 

determines the characteristic stochastic time Tstoch = 7t^ which now depends 
on the renormalized stochastic frequency z/j- = i/j'/2 + a'^/2vT and the renor- 
malized stochastic amplitude ay = ct/v^- 

4 Temperature-independence of rate constants 

Above examples show that the characteristic function F{t) exhibits an ex- 
ponential drop with the stochastic characteristic times r^tochi U = 1,2,...). 
The number and the form of the characteristic times are dictated by a precise 
stochastic process. [At a large difference between the amplitude and the fre- 
quency of the discrete process (DP or TP), the main drop is governed by the 
only decay parameters f l27|l or fl29|) .] Substituting the characteristic function 
F{t) in Eq. fl23|l and calculating the respective integral, one can achieve an 
analytic expression for the averaged transition rate constant. 

In the case of DP, the use of the exact expressions ( l25l) - (!27j) reduces the 
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averaged rate constant to (note 7 = 7^, v = vd) 

+ n{ux) L{'^,v]VLab + ■ (30) 
Here, the quantity x^ah — ^abd^^ appears as the phonon assisting inter-state 
couphng. In what follows, we rewrite the expression fl30l) in the form 

Kab = i \Vab\' r d^^-^{H^) + 1] ^(7, ^ab - CO) 



h? Jo h^u"^ 

+ n(w) L(7, z/; r^afe + } (31) 

where the coupling to a heat bath is reflected in the spectral function ( ITSl) 
whereas the influence of the stochastic field on the a — )■ 6 transition is concen- 
trated in the specific stochastic field generated (SFG) Lorentzian 

It is important to stress that the rate constant fl3Tl) includes a stochastic field 
(via the parameters u and 7 = a^/z/) in a non-perturbation manner. Therefore, 
it becomes possible to analyze various stochastic regimes of the transition 
processes in macromolecules dependending on the value of the generic field 
parameters u and 7. 

Before applying the theory to describe the conformational quasi-isoenergetic 
transitions, let us analyze a situation when \Qab ± wa| ^- Under this condi- 
tion, the SFG Lorentzian -^(7, u; Q) is reduced to a standard Lorentzian 

^(7'^) = -^T7^ (33) 

71 7^ -I- iZ^ 

with 7 = 7£), cf. Eq. ( 127|) . It is necessary to note that form (133|) is also valid 
in the case of TP, if only the function F{t) demonstrates a single exponential 
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drop with decay parameter 7 = 77-, cf. Eq. (129|) . In Lorenzian (133!) . quantity 
2/17 can be treated as the broadening of the molecular energy levels. If the 
parameter 7 is so small that Lorentzian (1331) appears as a narrow peak, one can 
calculate the rate constant taking a limit 7—7-0. This converts the Lorentzian 
into the Dirac's delta-function 6{fl) thus reducing the rate constant (|3T|) to 

Ka, = ^^4^ {[n{Qab) + imab) + r^(fi,„)^^(^],J } . (34) 

[In Eq. fl34|) . 6{Q) is the Heaviside unit function.] Note that in the case of DP, 
the form flM|) is realized only in the limit of high frequency stochastic field 
(when both inequalities, i^d ^ o'd and ^ \^ab ± i^a|, are satisfied for each 
bath mode A). In the case of TP, the same expression works under the condi- 
tion (z/y + (7^/i't) ^ \^ab ± wa|, i.e. in the two opposite limits of either high 
frequency stochastic field {ut ^ o't, the adiabatic regime) or low frequency 
stochastic field {ut ^ o't, the quasi-static regime). For exoergic transitions, 
when Qab > and exp (hQab/kBT) ^ 1, the rate constant ( l34l) is independent 
of temperature attaining the optical limit K^l^^^ = {\Vab\'^ / K^^ab) ^abi^ab) ■ 
Conversely, in the endoergic case when Qab < and exp {—h\Qab\/kBT) <^ 1, 
the rate constant appears as K^l^^^ = K^l^^^ exp {—E^'^^'^/kBT), thus demon- 
strating the Arrhenius like law with -E^b^*^ = ^|f^ab| being the activation energy. 

Now let us consider the formation of the rate constant at the quasi-isoenergetic 
transitions {Qab ~ 0). To this end, one has to analyze the action of the 
moderately high frequency stochastic field. For such field, the condition 7 ^ 
\Qab ± i^a| is satisfied, and thus, the SFG Lorentzian (|32i) is transformed into 

Lh,n^O) = —. (35) 
7r7 

Note a fundamental fact that at the quasi-isoenergetic transitions, forms ( l33l) 
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and (135|) are valid not only for the condition ufy ^ but for the condition 
a|) ^ as well. This means that during the DP, the expression (!35l) can 
be used to describe the quasi-iso energetic transitions at any relation between 
the amplitude and the frequency of the moderately high frequency stochastic 
field. Similar statement is true for the TP. 

The room temperatures correspond to the energies of the order 0.025 eV, iden- 
tical to 6-10^^s~^ or 200 cm~^. Therefore, if the transitions are accompanied by 
the vibrations of the order 60 cm~^ and lower, one can set n{ux) ~ ksT/hwx. 
This reduces the averaged rate constant fl3T]) to the form 



The physical origin of the stochastic parameter 7 is dictated by the random 
shifts of the molecular energy levels. We propose a phenomenological model 
where for the DP, the mean positive and the mean negative amplitudes (+o"d/2 
and — cr£)/2, respectively) are determined by the relation [{2n ha d)/ '2]'^ = SE'^ 
with SE"^ being the average of the square of energy fluctuations. In line with 
the general theory of thermodynamic fluctuations in canonic ensembles, this 
average is calculated as SE"^ = ksT'^idE / dT) [56] where E is the mean energy. 
It is possible to show that in thermodynamic equilibrium, the average linear 
frequency of the fluctuations is associated with the average energy E so that 
V£) = E/2Tih |15]. Thus, a stochastic parameter can be estimated with the 
relation 7^1 = {2k bT'^ / 7rh){d In E/dT). In the classical limit under consider- 
ation, the mean energy per a separate degree of freedom is equal exactly to 
thermal energy fc^T. Therefore, independently of the precise molecular sys- 
tem, the mean energy E is proportional to the bath temperature. It yields 
7£i = 2kBT/iTh. If TP is formed as the sum of the two DP, then Ut = 2z/£), 




(36) 
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O'T = ctd and thus ■Jt = (4/3)[(t|,^'d/(4z/|) + aj^)] = kBT/Snh. In summary, 
the fast equihbrium thermodynamic fluctuations form the stochastic decay 
parameter 7 = r^l^i^ that obeys a hnear dependence on temperature (in the 
classical limit). Thus, one can generally set 7 = C {2kBT /tiK) where factor 
is determined by the type of stochastic process. [In the above examples, one 
derives C = 1 and C = 1 /6 for the DP and the TP, respectively.] Consequently, 

Kab=^\Vab?Kh. (37) 

h 

In Eq. ( P7|) . Aafe = {2iiQ~^ duo {Jab{ou)/hPu^) is the temperature-independent 
factor. 

To demonstrate how the temperature dependence of the averaged rate con- 
stant (13T!) is canceled, we specify the form of a spectral function using the Song- 
Marcus's model [57f58] . For such model, Jab{^) = 27chw E^f^^"^^^ 6{u—uo) where 
ooq is the characteristic frequency of the bath phonon and Ej[f^^°^^^ = {gf^Y ^0 
is the reorganization energy of the a ^ b transition with participation of the 
only vibrational mode. The substitution Jahij^) in Eq. (13T|) reduces rate con- 
stant to the following reading form 

= iv^ + V+) ■ (38) 

Here, K^'' = {2TT\Vab\'^ /h^ool) E^fj'^"^^^ is the temperature-independent part of 
the rate whereas 

r]_ = u]o{l + n{uJo)) L{'y,u;nab-uJo) (39) 

and 

?7+ = Uq n{uJo) L(7, u; Vlab + ujq) (40) 
are the temperature-determining factors. 
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Fig. 1. Temperature behavior of the averaged transition rate constant (j38p at di- 
chotomous action of stochastic field. The physiologically important temperature 
region is presented. Curves 1 and 3 correspond to the off-resonant transition fre- 
quencies Qab = 10 cm~^ and flab = 80cm~^, respectively. Curve 2 manifests the 
cancellation of the temperature dependence at a resonant regime of transition at 
which transition frequency flab= 40 cm^^ coincides with phonon frequency ujq. The 
insertion shows the ratio r = rjj^/r]^ of the "off-resonant" and the "resonant" terms 
in the transition rate. The estimation of the factor rj = Kab/K^ = rj^ + r]^ is based 
on Eqs. 1^, ([Ml) and (00]) at 7 = 2kBT/-Kh and v = kBT/2T:h. 



In Fig. [H the temperature dependence of the transition rate constant is pre- 
sented. When the transition frequency VLab is far from the phonon frequency 
ujq ( "off-resonant" regime of transition) the factors ?7_ and 77+ are proportional 
to temperature (1 and 3 curves). The situation is changed significantly at the 
quasi-isoenergetic ("resonant") regime of transition when Vtab — wq- In this 
case, the main contribution to the rate constant is determined by the factor 
!]_ which now noticably exceeds the nonresonant factor 77+ (cf. insertion to 
Fig. [1]). But, at the quasi-isoenergetic transition, due to the linear tempera- 
ture rise of the stochastic parameter 7 the factor rj- demonstrates a practically 
temperature-independent behavior (compare curve 2 with curves 1 and 3). In 
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fact, one can talk about a temperature independent character of the transition 
process. 

5 Experimental support 

Eq. f l37|) and Fig. [1] present an important result, namely, that fast thermo- 
dynamic fluctuations are able to cancel the temperature dependence of rate 
constants characterizing the quasi-isoenergetic transitions in the macromolec- 
ular systems. In this section, we discuss the supporting experimental findings 
related to transfer processes in biological structures. In the begining of the Sec- 
tion 2, we have already noted the results concerning the observation of energy 
fluctuations in biomacromolecules. Here, we note an important finding that the 
fluorescence of the flavin cofactor indicates "all-or-none" enzymatic turnovers 
of the flavin reductase [38|39] while energy fluctuations in the pigment-protein 
comlexes of the light-harvesting antenna from photosynthetic bacteria exhibit 
the dichotomous behavior [ID|^ . It is quite possible to extrapolate these find- 
ings to other protein containing complexes. Thus, the kinetic equations f l22|) 
and respective rate constants fl3T|) derived by taking into account the averag- 
ing over the discrete energy fluctuations in a macromolecule, can be applied 
to explain the following experimental results. 

Example l-Temperature independence of desensitization onset of P 2X3 recep- 
tors. The P2X3 receptors belong to the family of ionotropic receptors widely 
evolved in the peripheral nervous system. These receptors are highly specific 
membrane proteins which link the binding of ATP molecules and their ana- 
logues to the opening and closing of the gate of a selective transmembrane 
ion pore [59]- The experimental data for the P2X3 receptors were obtained at 
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Fig. 2. Temperature independence of the desensitization onset of P2X3 receptors 
measured at physiological temperatures 25-40° C (selected from [19]) and its fit by 



phisiologically important temperatures of 25, 30, 35, and 40°C [19]. [In Fig|2l 
these data are reproduced with different symbols.] We have shown that the 
ATP-activated currents /(t)//(0) measured at different temperatures, mani- 
fest an identical two-stage decrease. To our knowledge, this is the first quan- 
titative observation of a temperature-independent gating in biological mem- 
branes. Evolution of the current is well described by the theoretical curve 
reflecting the double-exponential transition kinetics. 



with Pd{t) = 1 — /(t)//(0) being the desensitization probability of the chan- 
nels. In Eq. (141 p . Ai ^ 0.968 and A2 — 0.032 are the pre-exponential weights 
whereas ri ~ 14.7ms and T2 — 231ms are the temperature-independent char- 
acteristic times of the desensitization. 

Despite that primary, secondary and tertiary molecular structure of P2X3 
receptors is well-resolved [59], the particular molecular content of the recep- 
tor's gate(s) is yet unknown. But, it is possible to deduct that because of the 



Eq. gH). 



- A2 e-*/"2 , 



(41) 
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hydrogen and non-covalent (ionic) bond vibrations occur with characteristic 
frequencies 0.2-20 THz [60ll6T| . the fast local fluctuations of channel struc- 
tures can not directly lead to the slow open-closed transitions accompanied 
by extended flexible changes. However, these local fluctuations can create the 
stochastic fields varying the position of energy levels. 

We interpret the above experimental result in the framework of a model where 
the formation of an ion current appears as the sum of two partial currents. 
Let ij be the current through a separate open channel of the j{= 1, 2)th type. 
If Nj is the number of respective ion channels, then a partial current reads 
as Ij{t) = Njij Pj^^\t) where P^^\t) is the probability for the channel to be 
in an open state "o". Desensitization of each channel occurs independently 
reflecting the process of conformational transition from the open (conductive) 
state into the closed (nonconductive) state so that Pj^^\t) = exp(— t/rj). 
Introducing the quantities /(O) = Y,j=i.2 Ai = ^/(l + ^) and 

A2 = 1/(1 + .^) where C, = Niii /N2 12 , one expresses an ion current in the 
form I{t) = /(O) Poit) where Po(t) = Ej=i,2 Po \^) is the apparent (statis- 
tically averaged) probability of an ion channel to be in an open state. Bearing 
in mind that Pd{t) = 1 — Po{t), one arrives at the Eq. pil) (see also [62]). The 
physical explanation of the temperature independency of the desensitization 
process can be given based on the above proposed model of quasi-isoenergetic 
transitions in flexible molecules. Denoting via a = joa and b = jdjH the confor- 
mational isoenergetic substates participating in the open-close transition, we 
introduce the integral occupancies of states o and d as Po \t) = Y,a=i Pjoa{t) 
and Pjf\t) = I]/3=i Pjdpii) = 1 — Po\t) where /ijo and are the numbers 
of the respective degenerated substates for the jth type of channel. Eq. (l22l) 
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now reads 



(42) 



where 




(43) 



is the forward rate constant and (cf. Eq. (1371) ) 

Kj = — iK^r E E ^Joajdp . (44) 

''' a=l 13=1 

[The backward rate constant is given by expression K^^X^ = Kj/^jd] 

Since the disorder characteristic for desensitization increases the degener- 
acy of the molecular state, then fija ^ fijo- This reduces the kinetics of 
quasi-isoenergetic transitions to the nonrecurrent single-exponential kinet- 
ics with the characteristic time Tj ^ i^o-^d)"^ type of chan- 
nel. Thus, physically, the receptor desensitization appears as a temperature- 
independent open— )-close transition process between the degenerated quasi- 
isoenergetic molecular conformations of an ion channel. 

Example 2. -Temperature-independent degradation of an endogenous PER2. 
Very recently, a direct observation of temperature- independent conformational 
transformations in proteins has been provided in living clock cells [20]. It is 
known that a circadian periodicity in different organisms ranging from bac- 
teria to mammals remains robust over a wide range of temperatures (see e.g. 
|63]). One explanation of this characteristic may be that it is due to the tem- 
perature insensitivity of the degradation of PER2 protein which is the key 
period-determining protein in the mammalian circadian clock cascade. It has 
been shown that the phosphorylated PER2 degrades with a particular rate 
which was extremely sensitive to chemical perturbation, however, was remark- 
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Fig. 3. The temperature independence of the degradation of PER2 protein in mouse 
embrionic fibroblasts (reproduced from Figs. 2D and 2E of ref. [20J). The degrada- 
tion is monitored as a decay of luminescence (in percentages) of the PER2-luciferase 
complex. The diamonds indicate the data at 37° C and the line represents an ex- 
ponential curve. The insertion depicts the averaged life times calculated from the 
approximated curves at 27, 32 and 37° C. 

ably unaffected by physical perturbation such as temperature shifts ^U\. In 
vivo, the degradation of an endogenous PER2 is regulated by casein kinase 
le-dependent phosphorylation. During the phosphorylation, the ATP molecule 
is hydrolized thus donating inorganic phosphate oxygen groups to the chemi- 
cal structure of PER2. These groups, being negatively charged, can generate 
high frequency hydrogen bond fluctuations which, in turn, will affect the pro- 
tein degradation related to slow conformation motions. Therefore, as in the 
previous case of P2X3 receptors, the model of the stochastic regulation of the 
transition rate is also applicable. 



Experimentally, the PER2 degradation can be monitored by a decay of the bio- 
luminescence (cf. Fig. |3]) and appears as a single exponential drop of the ampli- 
tude of bioluminescence signal, A{t) = exp {—kdegt) where kdeg ~ 1-7 ■ lO^^s^^ 
is the slow degradation rate constant. As one can see, the half-times of the de- 



29 



cay curves has negligible temperature dependence (see the insertion to Figj3]) 
and may be considered temperature-insensitive. The physical mechanism of 
such temperature independency in physiologically important temperature re- 
gion, can be also associated with the equilibrium thermodynamic fluctuations 
and a pronounced disorder in the protein structure. This reduces the kinet- 
ics of the quasi-isoenergetic transitions to the nonrecurrent kinetics with the 
temperature-independent rate constant kdeg similar to that is given by Eq. 

m- 

6 Conclusions 

In this communication, a stochastic master equation (fTOl) for the diagonal part 
of the density matrix Pd(t) of an open quantum system is derived under the 
assumption that state energies Ea{t) are affected by random time-dependent 
fields. Special attention is given to derivation of kinetic equations for the aver- 
aged state occupancies of a macromolecule, Eqs. f|T^ and fl22|) . and respective 
averaged rate constants, Eq. fl2^ . The derivation is based on the model that 
incorporates two different types of motion in macromolecular system. The 
first type of motion is associated with a heat bath. The energy of the phonons 
belonging to a heat bath, covers the energy difference between the molec- 
ular states involved in the transition. Therefore, the bath phonons directly 
participate in the transition. The second one is characterized by the high fre- 
quency vibrations. The energy of the respective phonons strongly exceeds the 
transition energy, so that the high frequency vibrations do not directly par- 
ticipate in the transition. At the same time, the high frequency vibrations 
can create stochastic fields varying the energy of molecular states. The details 
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of the exact averaged procedure are demonstrated for two types of Kubo- 
Anderson stochastic processes, the dichotomous and the trichotomous ones. 
We found the characteristic times T^H^f^ of the drop of the averaged stochastic 
functional F{t). Additionally, the conditions when the F{t) exhibits a single- 
exponential drop, have been established. Special attention was focused on the 
analysis of the quasi-isoenergetic transitions. Here, a novel physical mecha- 
nism was proposed to explain the formation of the temperature-independent 
transition processes in molecular systems. In contrast with the well established 
mechanism describing a quantum phononless site-to-site low temperature par- 
ticle tunneling or a tunneling at room temperatures that is accompanied by 
an emission of high frequency phonons (when frequency Ux of the separate 
phonon satisfies the condition ux ^ ksT/h) piTU] ). the mechanism proposed 
here works in a classic (physiological) region of temperatures and at a weak 
coupling to the bath phonons. At such conditions, the transition is accom- 
panied by creation or annihilation of a single low frequency phonon of the 
frequency ux ^ ksT/h. The lack of temperature dependence in rate con- 
stant fl3Tl) occurs due to the thermodynamical stochastic variation of energy 
levels participating in the transitions. The frequency of this variation has to 
be of the order of the characteristic thermal frequency ksT/h. This reduces 
an averaged rate constant (I3T]) to a more simple form (136!) which, in turn, is 
simplified to the temperature-independent form (157|) . The two recent experi- 
mental observations of temperature-independent transitions in biosystems, are 
discussed. We show that these results can be explained in the framework of 
the model where the transitions are performed between the quasi-isoenergetic 
conformation states of macromolecules in the presence of random thermody- 
namic fluctuations of respective energy levels. In this case, both degradation 
and desensitization of proteins are generally accompanied by a pronounced 
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disorder in the protein structure thus essentially increasing the degeneracy 
of the quasi-isoenergetic molecular states. The explanation of the experimen- 
tal results is based on the model where part of molecular motions refer to 
fluctuations. Note that the presence of the stochastic motions in biological 
macromolecules is supported by spectroscopic measurements, thus indicating 
the existence of energy fluctuations of all-or-none type in protein containing 
complexes [38|39f40|HT| . 
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